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Under an in-plane magnetic field, the density of states of quasi two-dimensional carriers deviates 
from the occasionally stereotypic step-like form both quantitatively and qualitatively. For the first 
5— ( . time in the literature as far as we know, we study how this affects the spin-subband populations 

and the spin-polarization as functions of the temperature, T, and the in-plane magnetic field, B, 
[ for narrow to wide dilute-magnetic-semiconductor quantum wells. We examine a wide range of 

material and structural parameters, focusing on the quantum well width, the magnitude of the spin- 
spin exchange interaction, and the sheet carrier concentration. Generally, increasing T, the carrier 
spin-splitting, Uoa, decreases, augmenting the influence of the "minority" -spin carriers. Increasing 
, B, Uoa increases and accordingly carriers populate "majority" -spin subbands while they abandon 

^ ' "minority" -spin subbands. Furthermore, in line with the density of states modiflcation, all energeti- 

cally higher subbands become gradually depopulated. We also indicate the ranges where the system 
is completely spin-polarized. 
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I. PREAMBLE 



^ , An in-plane magnetic field applied to a quasi two-dimensional system distorts the equal-energy surfacesii^ or equiv- 
Q ' alently the density of states^i^ (DOS). An interplay between spatial and magnetic confinement is established and - 
I— ' ' properly- it is necessary to compute self-consistently the energy dispersion, -E^.o- (fc^; ) , where i is the subband index, a de- 
notes the spin, is the in-plane wave vector perpendicular to the in-plane magnetic field, B (applied along y), and z is 
the growth axis. Hence, the envelope function along z depends on i.e., il'i.a,k^,ky(j) oc Ci.<T,k^{z)exp{ikxx)exp{ikyy). 
This modification has been realized in magnetotransport^ and photoluminescence^ experiments. An impressive fluctu- 
00 ' ation of the in-plane magnetization in dilute-magnetic-semiconductor (DMS) structures in cases of strong competition 
. between spatial and magnetic confinement has been predicted at low enough temperatures^ and a compact DOS 
' formula holding for any type of interplay between spatial and magnetic confinement already exists^. 

Although this DOS modification can be extremely significant both quantitatively and qualitatively it is sometimes 
neglected without a second thought. Naturally, in the limit of very narrow quantum wells (QWs) or for B ^ 0, the 
DOS preserves the ideal step-like form. The "opposite" asymptotic limit is a simple saddle point, where the DOS 
"j^ diverges logarithmically'^. However, generally, the van Hove singularities which show up are not simple saddle points^. 

a Summarizing, models which ignore the above DOS modifications can only be applied to very narrow QWs or for 
, ; B^o. 

''O ' During the last years, the progress in growth, characterization and understanding of transition-metal-doped semicon- 
Ch , ductors has been impressive*'^' As a result, new phenomena have been discovered, e.g. tunnel magnetoresistance, 
■ spin-dependent scattering, interlayer coupling due to carrier polarization, electrical electron and hole spin injection, 
and electric field control of ferromagnetism^i^. Usually the host material is a III-V semiconducto r*:^: ^° . For example, 
in (Ga,Mn)As or in (In,Mn)As, Mn substitutes a small fraction of cations providing holes and local magnetic moments. 
Hence, the corresponding structures utilize the valence band. The highest ferromagnetic transition temperature, Tc, 
reported so far for Ill-V-based valence-band magnetic semiconductors is 173 K in (Ga,Mn)As epilayersi^. 
' In H-VI materials, Mn provides only local magnetic moments. The corresponding heterostructures, for example 
ZnSe/Zni_.a;_yCd2;MnySe, utilize either the conduction or the valence band, depending on the type of dopants used 
in the barriers, namely, donors (e.g. CI, I) or acceptors (e.g. Li), respectively. In the present article we investigate 
such a system where either the conduction band or the valence band can be exploited for spintronic applications. The 
key material of each structure (e.g. ZnSe, CdTe etc) may possess quite different material parameters e.g. positive 
or negative g factorsii. We also note that the band gap of common II- VI crystals covers all the range from the 
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ultraviolet to the infrared-^. Interestingly the existence of ferromagnetic order in n-doped (Cd,Mn)Te based structures 
-at extremely low temperatures- has been suggested both experimentally and theoreticall y "'^^i-'^^ . 

In a recent publication^ we restricted ourselves to DMS structures utilizing the conduction band and to very low 
temperatures. We studied the spin-subband populations, the internal and free energy, the Shannon entropy, and the 
in-plane magnetization M as functions of the in-plane magnetic field, for different degrees of spatial confinement. The 
enhanced electron spin-splitting Uoa can be considered as the sum of two terms, a and [3. a is proportional to the 
cyclotron gap, hwc, while (3 arises from the exchange interaction between the itinerant carrier (conduction electron 
in ReL^) and the localized spins (Mn+^ cations in RefJ^). Notice that in such an approximation the direct exchange 
interaction between the neighboring localized impurity spins is neglected, being much smaller than the interaction 
between impurity spins and carrier spinsi^, although according to a recent report it might influence the carrier spin 
polarization^*^. The very low T impelled us to a drastic first approximation^, i.e. to take into account only (3, and 
moreover to approximate the corresponding Brillouin function by 1. 

In the present article we attempt a major improvement. Namely, we examine the relative influence of a and /3 in a 
wide temperature band (0 to 400 K) and in a wide in-plane magnetic field band (0 to 20 T), as well as in a wide range of 
material parameters, not necessarily restricting ourselves in the conduction band^^. Our purpose is to systematically 
study the influence of the DOS modification on the spin-subband populations and the spin-polarization of quasi two- 
dimensional carriers, as functions of the in-plane magnetic field and the temperature. Besides, we indicate the ranges 
where the system is completely spin-polarized. In Section|IT]we introduce our theoretical framework. In Section Hill we 
examine the spin-subband populations and the spin polarization, C, of non-magnetic-semiconductor (NMS) / narrow 
to wide dilute-magnetic-semiconductor (DMS) / NMS quantum wells (QWs), as a function of the temperature, T, and 
the in-plane magnetic field, B. We notice that in the present system due to the influence of carriers, increase of the 
QW width transforms the heterostructure from an "almost perfect square QW" to a "double QW with a soft barrier" 
("a system of two separated heterojunctions")^. Thus, the present heterostructure allows us to study "single" as well 
as "double" QWs. To facilitate the reader, we provide in Fig. [1] sketches of the self-consistent QW profiles for T = 
20 K, i? = T and sheet carrier concentration Ns ~ 1.566 x 10^^ cm"^, for QW widths 10 nm, 30 nm and 60 nm. 
We examine how the DOS modification affects ( for a wide range of material and structural parameters focusing on 
the quantum well width, the magnitude of the spin-spin exchange interaction coupling strength, and the sheet carrier 
concentration. Finally, in Section [TVl we briefly state our conclusions. 
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FIG. 1: (Color online) Sketch of the QW profiles for T = 20 K, B = T and A^^, = 1.566 x lO" cm-^ for QW widths 10 nm, 
30 nm and 60 nm. 



II. THEORY 



Under a magnetic field, B, applied parallel to the interfaces, the equal energy surfaces are gradually distorted. The 
density of states deviates from the ideal step-like form both quantitatively and qualitatively^, i.e. it takes the form: 



(1) 
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where it is implied that the QW is along the z axis and the magnetic field is applied along the y axis. 8 is the 
step function, A is the xy area of the structure, to* is the effective mass^l. Ei^a^^x) are the spin-dependent xz-plane 
eigenenergies. Generally, Ei^a{kx) must be self-consistently calculatedi^i^ii^i^i^ii. Equation H]) is valid for any type 
of interplay between spatial and magnetic confinement. The dependence in Eq. ([1]) increases the numerical cost 
by a factor of 10^ — 10^ in many cases. This dependence is quite often "conveniently" ignored, although this is 
only justified for narrow QWs. However, with the existing computing power, such a "simplification" is not any more 
necessary. Only in the limit i? — > 0, the DOS retains the occasionally stereotypic staircase shape with the ideal step 
^■^^-T- for each spin. The opposite asymptotic limit of Eq. ([T|) is that of a simple saddle point, where the DOS diverges 
logarithmically^. The DOS modification significantly affects the physical propertiesii^ii^i^i^i^ii. For completeness, we 
notice that Eq. ([T]) ignores the effect of disorder which -with the current epitaxial techniques- is important when the 
concentration of magnetic ions is higb^i^. Disorder will certainly induce some broadening of the spin-subbands. 

In DMS structures, the electron spin-splitting, Uoai is not proportional to the cyclotron gap, hiUc, i.e. it acquires 
the form2^^^: 



Uoa = ^i^f^c - yNoJsp-dSBsiO ^a + p. (2) 

a = a{B) describes the Zeeman coupling between the spin of the itinerant carrier and the magnetic field, while 
/3 = P{B,T) expresses the exchange interaction between the spins of the Mn+^ cations and the spin of the itinerant 
carrier (initially supposed to be an electron), g* is the g-factor^^ of the itinerant carrier, y is the molecular fraction 
of Mn. A'o is the three-dimensional (volume) concentration of cations. Jsp~d is the coupling strength due to the 
spin-spin exchange interaction between the d electrons of the Mn+^ cations and the s- or p-band electrons, and it is 
negative for conduction band electrons. The factor SBs{C) represents the spin polarization of the Mn+^ cations. The 
spin of the Mn+^ cation is 5 = 5/2. i?s(C) is the standard Brillouin function, whil o^^'^^ 



gMntisSB — Jsp-dS- 



, 

ks is the Boltzmann constant, /i^ is the Bohr magneton. gMn is the g factor of Mn'^^. nciawn and n„p are the 
spin-down and spin-up three-dimensional (volume) concentrations measured e.g. in cm^'^, while Ns^down and N. 



s,up 



used below are the spin-down and spin-up two-dimensional (sheet) concentrations measured e.g. in cm~^. In Eq. [3] 
(and only there) we approximate ridown — riup ~ {Ns,down — Ns^up)/L, where L is the QW width. The first term in the 
numerator of Eq. [3] represents the contribution of the Zeeman coupling between the localized spin and the magnetic 
field. The second term in the numerator of Eq. [3] (sometimes called "feedback mechanism") represents the kinetic 
exchange contribution which -in principle- can induce spontaneous spin-polarization i.e. in the absence of an external 
magnetic field^S. Notice that Udown — nup is positive for conduction band electrons. Finally, for conduction band 
electrons, the spin polarization can be defined by 



-^s.down ^s,up / .\ 



Ns = Ng ,io^n + ^s.up is the free carrier two-dimensional (sheet) concentration. 

The use of such a simplified Brillouin-function approach is quite common when dealing with quasi two-dimensional 
systemsi2iii^i22ii2Si^. This way, the spin-orbit coupling is not taken into account. This is certainly a simplification, 
since increasing temperature, the magnetization of the magnetic ions competes with spin-orbit coupling. The spin- 
orbit coupling^^'^^ induces temperature dependent spin relaxation. Therefore, the carriers' spin-polarization docs not 
only depend on the magnetic order of the magnetic ions, expressed here with the help of the Brillouin function and 
the carriers' spin relaxation influences the magnetic order of the localized magnetic moments. 

The variation of the temperature, T, affects the spin polarization. The spin polarization is also influenced by the 
magnetic fleld, in an opposite manner i.e. B tends to align the spins. Furthermore, for each type of spin population, the 
in-plane magnetic fleld -via the distortion of the DOS- redistributes the electrons between the subbands. Consequently, 
the spin polarization can be tuned by varying the temperature and the magnetic field. Indeed, preliminary conduction 
band calculations for specific values of the material parameters, for very narrow quantum wells, have shown^ that 
when the "feedback mechanism" due to the difference between the populations of the spin down and the spin up 
electrons can be neglected, the spin polarization vanishes for _B — > 0. The analysis presented above can be useful for 
p-doped structures, assuming -as usual- that a single valence band description is a fair first approximation. 
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III. RESULTS AND DISCUSSION 



Initially we consider heterostructures of the type n-doped ZnSe / Zni_2;_j,Cda;Mnj^Se /n-doped ZnSe. Let us take 
y = 0.035, —yNQJsp_d ~ 0.13 Hartree*, and the conduction band offset, AUct = 1 Hartree*— . We notice that for 
ZnSe, 1 Hartree* « 70.5 meV. ZnSe has a sphalerite-type structure and the lattice constant is ~ 0.567 nm. Hence, 
— Jsp-d ~ 12 X 10~^ eV nm'^. This is one order of magnitude smaller than the value commonly used for the III-V 
Ga(Mn)As valence band system {Jp^ 



= 15 X 10-2 eV nm3)ia^>^. 
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FIG. 2: (Color online) The spin-subband populations Nij = Nij(B) for T = 20 K (first row), and Ntj = Nij{T) for B = 
10 T (second row) of a n-doped ZnSe / Zni-^-yCd^MnySe /n-doped ZnSe QW with y = 0.035. -Jsp-d = 12 x 10"^ eV 
nm"^. 00 stands for the ground-state spin-down-subband, 10 for the 1st excited spin-down-subband, 01 for the ground-state 
spin-up-subband, and 11 represents the 1st excited spin-up-subband. Each column corresponds to different well width i.e. L = 
10 nm, 30 nm, and 60 nm. 

Figure depicts the spin-subband populations, Nij as a function of B (a-c), and as a function of T (d-f), for 
three different well widths, namely (a,d) for .L = 10 nm (b,e) for L = 30 nm, and (c,f) for L — 60 nm. Initially, we 
deliberately keep the total sheet carrier concentration constant {Ng ~ 1.566 x 10"'^^ cm"^), assuming that all dopants 
are ionized. In (a-c) T = 20 K. In (d-f) .B = 10 T. The pair ij is defined in the following manner: 00 symbolizes the 
ground-state spin-down-subband, 10 the 1st excited spin-down-subband, 01 the ground-state spin-up-subband, and 
finally 11 symbolizes the 1st excited spin-up-subband. Due to the small value of Jsp-d, the influence of the "feed-back 
mechanism" due to the difference between spin-down and spin-up concentrations is negligible in the present system. 
Indeed, since — J.sp^d "'^"""^"""'' negligible here, then for B = 0, it follows that (a) C ~ 0, thus -Bs(C) ~ 0: therefore 
/3 « 0, and (b) g*fiBB — a — 0. Hence, Uoa ~ 0, and consequently C ~ 0. In fact, inspection of Figsl2] (a-c), 
reveals that for S = 0, in FigEfa) iVoo = ^^oi, in FigH^b) iVoo = ^oi and = Nn, and in Figl^l^c) iVoo = ^oi 
and A'^io — Nn. For the very wide quantum well {L = 60 nm), as expected^"^, the four spin-subbands are almost 
equally populated for .6 = 0. Increasing B, we observe that there are two mechanisms which cause depopulations: 
(I) The increase of Uoa- eliminates spin-up electrons, namely A^oi and A^n continuously decrease, increasing B. (II) 
The DOS modification which depopulates all excited states, regardless of their spin"''^, namely the eventual decay of 
iViQ. Finally, in FigsH] (d-f), we witness the survival of only A^oo at very low T, since Uoa acquires its bigger value at 
zero temperature. Increasing T, Uoa- decreases, augmenting the influence of the spin-up electrons. 

Figure [3] depicts the relative influence of the Zeeman term, a, and the exchange term, /3, in wide B and T ranges, 
for a n-doped ZnSe / Zni_a;_j,Cd2:MnySe /n-doped ZnSe QW with .L = 60 nm and y — 0.035. In reality L is of 
no importance here due to the negligible impact of the "feedback mechanism" with these material parameters. For 
comparison we notice that the conduction band-offset, AUcb — 1 Hartree* « 70.5 meV. The spin splitting in the 
present article, Uoa = a + P, while U®^ was used in our previous low-T calculations^ (-^5/2(0 approximated by 1, and 
a ignored). Figure [3] elaborates the competition between B (aligning spins) and T (bringing on anarchy). Figure 
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FIG. 3: (Color online) The relative influence of the Zeeman term, a, and the exchange term l3 in wide B and T ranges, for a 
n-doped ZnSe / Zni_^_„Cd^MnySe /n-doped ZnSe QW with L = 60 nm and y = 0.035. -J^p-d = 12 x 10"^ eV nm^. Of 
course, a = a{B), while /? = f3{B,T). (a) a = a{B), P = [3(B) for T = 20 K. (b) a = a{T) = constant, (3 = /3(T) for B = 10 
T. L = 60 nm. Each panel also contains the spin-splitting, Uoa = Q + /3, as well as the value of of the spin-splitting used in 
our previous low-T calculation- (i.e. taking into account only (3 and approximating the corresponding Brillouin function by 1), 
U%. For comparison we notice that the conduction band-offset, Af/ct = 1 Hartree* ~ 70.5 meV. 




FIG. 4: (Color online) The spin polarization, C,, tuned by varying: (a) the in-plane magnetic field, B, keeping T = 20 K (left 
panel), and (b) the temperature, T, keeping _B = 10 T (right panel), for different well widths, L = 10 nm, 30 nm, and 60 nm. 
-Jap-d = 12 X 10"^ eV nm^ 



justifies our previous low-temperature approximation: at low enough T, Uoa ~ C^Jo-. At higher temperatures, B<^/2{0 
cannot be approximated with 1. As ksT increases, ^ decreases, and consequently i?5/2(0 < 1- In other words, 
increasing T, the spin-splitting decreases allowing enhanced contribution of the spin-up electrons to the system's 
properties. Finally we notice that an opposite sign of g* (e.g. CdTe vs. ZnSe) is expected to have small effect on the 
results since the most important term is [3. 

Figure m depicts the spin polarization tuned by varying the parallel magnetic field and the temperature, for different 
choices of the well width. Since for i? > 8 T, ^ = 1, only the range Be [0, 8 T] is presented in Fig. [3^. Since for T > 
150 K, C is less than « 0.1, only the range T £ [0, 150 K] is presented in Fig.[4j3. Because of the DOS modification^, 
resulting in different distribution of electrons among the spin-subbands (cf. Fig[2]), we observe a clear dependence 
of C = C{L), i.e. ({L — 60 nm) > ({L = 30 nm) > ({L = 10 nm). We also observe that for _B = 0, C vanishes, 
i.e. there is no spontaneous spin polarization phase due to the tiny "feedback mechanism" for this choice of material 
parameters. 

Subsequently we deliberately increase —Jsp-d by an order of magnitude i.e. we present in Fig. [5] results with 
— Jsp-d = 12 X 10~^ eV nm'^ (which is of a little smaller magnitude than the value commonly usedi^i^^!^ for the 
III-V Ga(Mn)As valence band system, Jpd = 15 x 10~^ eV nm'^). L = 60 nm and T = 20 K. Comparing Fig. [5] with 
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FIG. 5: (Color online) The spin-subband populations, Nij (left panel) and the spin polarization, ( (right panel) tuned by 
varying B for L = 60 nm, T = 20 K, using -Jsp-d = 12 x 10"^ eV nm^. 




FIG. 6: (Color online) The spin-subband populations, Nij (left panel) and the spin polarization, (right panel), tuned by 
varying the sheet carrier concentration, Ns, for L = 60 nm, T = 20 K and B = 0.01 T, using J = 12 x 10~^ eV nm"^. The little 
arrows indicate Ns values where we also compare with B = 0.0001 T in the text. 



Fig. [2t and Fig. 2^ we observe that: {a') The greater value of —Jsp^d makes it much easier to attain a completely 
spin-polarized system (C = 1) i.e. for i? > 1 T instead of B > 8 T. (/?') Initially, increasing B, due to the increased 
Uoa, NiQ grows, in contrast to Fig. Naturally, subsequently A^io is depopulated because of the in-plane magnetic 
field induced DOS modification. (7') Although the system is more susceptible to spin-polarization, still, practically 
no spontaneous spin-polarization phase exists for i? = 0, at this temperature. 

Up to now, we have deliberately kept the total sheet carrier concentration constant. Below we examine the influence 
of Ns on the spin-subband populations and the spin polarization for different values of the magnitude of the spin-spin 
exchange interaction, J. Since Ns is affected by many factors (QW profile, material properties, valence-band- or 
conduction-band-based structures etc) we have decided to use J as a parameter here. Naturally, in a heterostructure 
where higher Ng can be achieved we may require smaller values of J in order to completely spin-polarize carriers. 
Using the rest material parameters as above but modifying J, we have systematically studied the Ns influence. For 
J = 12 X 10~^ eV nm'^ there is a very small influence of Ns on C. The situation changes using J = 12 x 10^^ cV 
nm^. Figure [H shows Nij and C tuned by varying Ns for L = 60 nm, T = 20 K and B — 0.01 T, using J = 12 x 10~^ 
eV nm'^. We observe that increase of Ns from « 1.0 x 10^ cm~^ to « 1.0 x 10^^ cm^^ is sufficient to completely 
spin-polarize carriers. This is purely due to the "feedback mechanism" stemming from the difference between the 
populations of spin-down and spin-up carriers. If we decrease B from 0.01 T to 0.0001 T, then e.g. (a) for Ns = 1.175 
X 10^ cm-2, ( changes from 0.497 to 0.005, (b) for Ns = 3.917 x 10^° cm^^ ( changes from 0.973 to 0.909, however, 
(c) for Ns = 1.175 X 10^^ cm~^, ( remains 1. 
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IV. CONCLUSION 



We have studied the spin-subband structure of quasi two-dimensional carriers in dilute-magnetic-semiconductor- 
based heterostructures, under the influence of an in-plane magnetic field. The proper density of states was used for 
the first time, incorporating the dependence on the in-plane wave vector perpendicular to the in-plane magnetic 
field. We have examined the interplay between different degrees of spatial and magnetic confinement, as well as the 
influence of temperature in a wide range. We have systematically studied the spin-subband populations and the 
spin-polarization as functions of the temperature and the in-plane magnetic field. We have examined a wide range of 
material and structural parameters, focusing on the quantum well width, the magnitude of the spin-spin exchange 
interaction, and the sheet carrier concentration. In particular we have shown that with sufficient magnitude of the 
spin-spin exchange interaction, the sheet carrier concentration emerges as an important factor to manipulate the 
spin-polarization, inducing spontaneous (i.e. for B ^ 0) spin-polarization. We have shown how at low temperatures 
the spin-splitting acquires its bigger value and how it decreases at higher temperatures. Increasing the in-plane 
magnetic field, the spin-splitting increases inducing depopulations of the "minority" -spin subbands. Moreover, the 
DOS modification induces depopulations of all energetically higher subbands. Finally, we have indicated the ranges 
where the system is completely spin-polarized. 
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